超 临 齐 水 传 热 关联 式 预测 传 热 恶化 的 评价 
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摘 要 : 本 文 在 较 宽广 的 工 况 范围 内 收集 、 整 理 了 垂直 上 升 光 管内 超 临 界 压力 水 的 对 流传 热 实验 数据 (包括 传 热 恶 
化 实验 数据 ) 以 及 超 临界 压力 水 传 热 的 计算 关联 式 。 将 收集 到 的 传 热 关联 式 与 传 热 恶化 实验 数据 进行 比较 ， 进 行 误 
差分 析 。 对 传 热 关联 式 的 形式 分 析 后 ， 对 垂直 上 升 管 中 传 热 实验 数据 进行 拟 合 得 到 一 个 新 的 传 热 关 联 式 并 对 其 进行 
评价 ， 效 果 较 好 。 
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Evaluation of supercritical-water heat-transfer correlations for 


heat transfer deterioration 
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Abstract: The available experimental data of heat transfer for water and the empirical correlations which 


are used to predict heat transfer at supercritical pressure in vertical upward tubes are extensively collected. 
The calculated values of these experimental correlations are compared with the experimental value of the 
heat transfer deterioration, errors and accuracy analyses for these correlations arrived are presented. The 
heat transfer experimental data in vertical upwards tubes are fitted to get a new heat transfer correlation. 
An error analysis for the new correlation is presented which works well with heat transfer deterioration. 
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0 引言 


超 临 界 水 冷 堆 《SCWR) 作为 第 四 代 核 能 系统 超 临界 压力 水 的 对 流传 热 实验 数据 《包括 恶化 实验 
中 唯一 的 水 冷 堆 ， 结 合 了 世界 上 大 量 应 用 的 轻 水 堆 数据 ) 以 及 预测 管内 的 超 临界 压力 水 传 热 的 计算 关 
以 及 超 临 界 火电 技术 的 特点 ,热效率 达到 40% 以 上 ， 联 式 。 本 文 在 收集 的 传 热 恶 化 实验 数据 基础 上 ， 评 
其 系统 比较 简单 ， 是 最 有 可 能 成 为 以 后 生产 低 成 本 价 了 传 热 关 联 式 对 传 热 恶化 的 预测 能 力 。 通 过 对 超 
电力 的 主力 堆 型 门 。 管 内 的 超 临 界 压力 水 传 热 特性 临界 压力 水 传 热 关 联 式 的 分 析 ， 对 收集 的 传 热 实验 


的 研究 是 研发 超 临界 机 组 和 超 临 界 水 冷 堆 的 基础 ， ”数据 进行 线性 回归 ， 得 到 一 个 新 的 关联 式 ， 其 对 超 
把 超 临界 水 作为 工作 流体 最 关心 的 就 是 它 在 拟 临 界 。 ”临界 水 的 传 热 恶 化 现象 有 一 定 的 预测 能 
点 附近 物性 的 剧烈 变化 。 
热 实验 数据 信息 以 及 传 执 
现 有 实验 研究 发 现 ， 根 据 热流 密度 和 质量 流速 | 伟 热 实验 数据 信息 以 及 传 热 天 联 式 
的 不 同 出 现 了 3 种 类 型 的 传 热 现象 四。 传 热 正常 通过 文献 调研 ， 广 泛 收集 整理 了 超 临界 压力 重 


主要 发 生 在 远离 临界 或 拟 临界 点 区 域 ， 其 传 热 系数 。 直上 升 光 管 传 热 实 验 工 况 202 个 ， 共 10655 个 数据 


可 用 Dittus-Boelter 等 一 类 传 热 关联 式 较 好 预测 ; 传 点 B21， 其 中 包括 了 34 组 传 热 恶化 实验 工 况 ，1164 
热 强化 ， 发 生 传 热 强化 时 ， 传 热 系数 在 拟 临 界 点 附 。 人 他 热 恶化 实验 数据 点 

近 有 传 热 系数 的 局 部 峰值 ， 传 热 恶化 ， 传 热 系数 在 NUR xL. e USE 
BUSZLRHDERGIGRMUH RARA BARERA LL 二 个 和 "oHm 
升 高 。 由 于 水 物性 在 拟 临 界 点 附近 的 剧烈 变化 ， 超 Lo VURU DERE TASEBTR UR 7] P RURGRELG 
临界 压力 水 传 热情 况 就 变 得 非常 复杂 ， 传 热 恶化 现 。 关系 式 。 其 中 绝 大 多 数 都 是 建立 在 变 物性 的 单 相 强 


象 就 需要 被 细致 的 研究 。 制 对 流传 热机 理 上 ， 大 部 分 采用 努 谢 尔 特 数 、 雷 诺 
本 文 收集 、 整 理 了 公开 文献 中 垂直 上 升 光 管内 数 、 普 朗 特 数 及 其 它 物性 参数 比值 形式 的 修正 项 诸 


201802.00076v1 


chinaXiv 


Ülu/uo k/k KWE KRA FE EE 
般 采 用 流体 温度 、 壁 面 温度 或 膜 温度 中 的 一 种 。 本 
文 在 公开 文献 中 ,收集 了 21 个 垂直 上 升 管内 超 临界 
压力 水 传 热 的 计算 关联 式 呈 9。 


2 不 同 关 联 式 预测 传 热 恶 化 能 力 评价 


图 1 给 出 了 不 同 关联 式 传 热 系 数 的 计算 值 与 
胡 志 宏 中 "实验 值 的 比较 结果 。 实 验 参 数 是 ，G=600 
kg/m’s, q=300 kW/m?, P=30 MPa, D=26 mm. Dittus 
and Boelter(1930)5 关 联 式 在 传 热 恶化 区 域 , 传 热 系 
数 先 上 升 到 一 个 局 部 峰值 后 再 下 降 ， 在 传 热 恶化 区 
域 与 实验 结果 相反 。 其 余 关 联 式 在 一 定 程度 上 预测 
到 了 传 热 恶化 现象 ， 但 是 预测 的 结果 差异 较 大 。 

为 了 进一步 评价 不 同 关联 式 ， 通 过 平均 相对 误 
28 0 人、 平均 相对 绝对 值 误差 co* 和 标准 偏差 0 ,三 个 
统计 指标 定量 评价 关联 式 的 精确 度 。 具 体 定义 如 下 : 


0, -yey/N (1) 
i=] 
c, - V. |e |/N (2) 
ial 
o, = Ye, -0,7 /(N -1) (3) 


i=l 
其 中 N 是 数据 点 个 数 , 相 对 误差 e; 由 下 式 计 算 : 

e, 2 [(Nu), — (Nu), ]x100 / (Nu),, (4) 
0, 表示 经 验 关 联 式 预测 值 和 实验 值 之 间 误 差 的 
总 体 趋势 ，a 消除 a 中 正 负 抵消 的 影响 ，cr, 表示 

预测 值 和 实验 值 的 分 散 度 。 
表 1 给 出 了 不 同 关联 式 预 测 值 与 传 热 恶 化 实验 
的 误差 分 析 。 不 同 指标 下 ， 每 个 指标 对 应 的 较 好 关 
联 式 并 不 一 致 。 表 3 对 三 个 指标 都 较 好 的 6 个 关联 式 
都 进行 了 加 粗 标识 。 这 5 个 关联 式 是 : Watts and Chou 
(1982)P"., Bishop(1964)7?l, 48i% (2004)P%, Jiyang 
Yu(2009)P8.. Mokry(201 1)? 5130] (2012)P " A BEC 
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1 关联 式 计算 传 热 系数 与 胡 志 宏 实验 值 "的 比较 


Fig 1 Comparison of heat transfer correlations with 
[17] 
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experimental data (Hu Zhihong 


3 新 的 关联 式 建立 
3.1 计算 关联 式 形式 讨论 

由 于 工 质 物 性 随 温 度 的 变化 ， 既 影响 流体 在 热 
流 方向 上 的 动量 交换 与 能 量 交 换 ， 也 影响 流体 在 热 
流 方 向 的 重力 场 ， 前 者 被 视 为 “ 变 物性 ”的 影响 ， 
后 者 称 之 为 “浮力 ”的 影响 的 。X.Chengi 总 结 了 超 
临界 压力 下 流体 传 热 关 联 式 ， 考 虑 流体 在 拟 临 界 附 
近 区 域 物性 剧烈 变化 对 换 热 造成 的 影响 ， 大 多 数 经 
验 关联 式 都 是 以 Dittus-Boelter 关 联 式 为 基础 进行 修 
正 得 到 。 


case) 


zr 


ZF 


Nu, = Cx Re?x Pr" x F (5) 
修正 因子 Z 常 采用 壁面 温度 和 主流 体温 度 下 物 

性 比 形式 的 修正 方法 。 
Fa| p/P A Ag] (6) 


表 1 不 同 关 联 式 和 传 热 恶化 实验 数据 比较 的 误差 分 析 
Table 1 Comparison of data fit with deteriorated heat transfer 


experimental data 


作者 关联 式 u 95 E 


Dittus and Boelter1930??! — 187.449 189.908 — 281.696 
Bringer and Smith 19577?) 87.615 — 93.106 — 105.159 
Miropol'skiy & Shitsman 

91.646 94.615 107274 
1957P% 
Bishop 19642”! 24.254 28.752 39212 
Swenson 1965] 207.897 208.650 154.880 


Krasnoschekov 19672” 5.447 38.572 48011 


Ornatsky 19707?! 35.005 | 46.700 67211 
Yamagata 1972P” 52.036 55.427 59277 
Jackson 1975P?! 45.867 48.968 55.971 
Jackson and Fewster 

197502 39259 41.974 49.457 
Watts and Chou 19827" 20.890 28.419 39.494 
Razummovskiy 199021 116.102 118.584 — 119.443 
Gorban 19990?! 26.991 43446 . 79.459 
TAE NUERA 20.754 25.963 40.197 
定义 书签 。] 

徐 峰 20042" -16.093 — 30.293 32.764 
Bo 2008P?! 19.999 34.180 — 47.994 
Jiyang Yu 200929 23.398 — 28.040 37.958 
朱 晓 静 2009P" -29.384 40.059 37210 
Mokry 2011/7! 0.041 19.04 30.190 
Xg 2012/7! 4.364 18.508 28.982 
FRA 291207! 60.892 65.488 58.401 


管子 在 加 热 过 程 中 ， 温 度 差 产 生 的 密度 差 会 导 


致 浮力 作用 影响 超 临 界 水 在 径 对 
Red d ^h. apis 
一 定 的 条 
^ mE EJES 
布 进而 改变 能 量 扩散 强度 ， 也 会 对 换 热 造成 一 定 的 
影响 。 文 献 记 -9 对 浮 升 力 和 热 加 速 对 换 热 的 影响 进 
行 了 定量 的 分 析 ， 用 6G，g 分 布 反映 浮力 作用 和 加 
速 作 用 对 换 热 模型 的 影响 ， 参 数 定义 为 : 


. Pgd'q, . dq, 
G 2 —————w WW. 


Av? GC 
因此 最 终 的 关联 式 形式 : 


zc A ý P 和 *\M x 
Nu=CRe" Pro| —| |=] (G) q™ (9) 
i (4 | (2 | ( ) 


3. 2 新 的 关联 式 拟 合 及 对 传 热 恶 化 预测 能 力 的 评价 


pb 


为 了 能 够 得 到 适用 范围 较 广 的 换 热 关联 式 ， 以 
整个 超 临界 垂直 上 升 光 管 传 热 实验 数据 源 〈 共 
10655 个 数据 点 ， 包 含 1164 个 传 热 恶化 数据 点 ) 为 
基础 采用 多 元 线性 回归 得 到 的 新 的 关联 式 如 下 : 


Nu = 0.032089Re088832 Pr; 71399 (4 nes 
b Ww 
x(p, / p, aid [on t quem T 
其 中 Pr = C,H c Gl, BH 
ZN , b 一 , — 
A * (T,-L) 


本 文 的 关联 式 平均 相对 误差 c |、 平均 相对 绝对 
值 误 差 o。 和 标准 偏差 o 3 分 别 是 -6.103，18.063， 
21.385， 标 准 偏差 最 小 ， 结 果 分 散 度 较 小 。 

本 文 关 联 式 与 其 他 5 个 关联 式 的 预测 值 与 传 热 
恶化 实验 值 对 比 结果 如 图 2 所 示 。Watts and 
Chou(1982)E9、Jiyang Yu(2009)P9、Mokry(201103] 
以 及 刘 奢 (2012)P 关 联 式 整体 预测 值 偏 高 ， 徐 峰 
(2004) POK SX Ty UH YE E + 309 EDGE V 2E S ER P 
较 多 ; 本 文 关 联 式 预测 值 分 散 度 最 小 。 

表 2 给 出 了 5 个 预测 效果 较 好 的 关联 式 与 本 文 
关联 式 计算 的 预测 值 与 实验 值 误差 分 布 在 10%， 
20% 和 30% 三 个 误差 带 内 的 数据 百分比 。 可 以 看 出 
Mokry(2011)?., j£ (2012) P p t A Sc BO SERES 
预测 值 70% 以 上 数据 落 在 土 30% 范 围 内 。 通 过 上 面 
的 分 析 ， 可 以 看 出 MokryQ011)?], xiii (2012) P 
以 及 本 文 的 关联 式 对 于 传 热 恶 化 的 预测 效果 较 好 。 
表格 2 不 同 误差 范围 关联 式 预测 的 数据 量 及 数据 量 百分比 


Table 2 Percentage of calculated-experimented Nu that drops in 


certain error bands 


作者 关联 式 10%〔 个 数 / 20% (个 数 / 30% 《个 数 / 
百分比 》 百分比 ) 百分比 ) 
Bishopl964U7 327/28.09%  617/53.0196 810/69.59% 
Watts and 
hp 304/26.12% 617/53.01% 843/72.42% 
徐 峰 2004P?! 196/16.84% — 404/34.7196 — 664/57.0494 
Jiyang 2009P9! — 323/27.7594  638/54.8199 — 819/70.3696 
Mokry 
2913 458/39.359, 791/67.96% — 992/85.2294 
XÆ 20127" ^ 469/40.299,  819/70.3696 1031/88.5796 
本 文 关联 式 370/31.79% 725/62.29% 918/78.87% 
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图 2 不 同 关联 式 预 测 值 与 传 热 恶 化 实验 数据 的 比较 
Fig 2 Comparison of experimental data (heat transfer deterioration) with predictions by different correlations 
4 LES [3] Griem H. A New Procedure for the Prediction of Forced 
ZA V 


本 文 在 收集 的 垂直 上 升 管内 超 临 界 水 换 热 实验 
数据 基础 上 , 评价 了 收集 的 21 个 换 热 关联 式 并 在 数 
据 基础 上 拟 合 了 新 的 关联 式 ， 得 出 以 下 主要 结论 : 

各 个 关联 式 的 预测 传 热 恶化 效果 有 很 大 差异 ,5 
个 有 一 定 的 预测 传 热 恶 化 能 力 的 关联 式 : Watts and 
Chou(1982) 、 徐 峰 (2004) 、 Jiyang Yu(2009) ~ 
Mokry(2011) 0 & x $8 (2012) IG 

本 文 考虑 了 物性 参数 、 浮 升力 以 及 热 加 速 对 传 
热 的 影响 ， 在 收集 到 传 热 实验 数据 基础 上 ， 回 归 得 
到 一 个 新 的 无 量 纲 的 传 热 关联 式 。 该 关联 式 与 
Mokry(2011) fl xi E (2012) 2& WX XU 4 Tf B5] RU s 
热 恶 化 能 

超 临界 换 热 关联 式 有 很 大 的 应 用 局 限 性 , 不 同 
关联 式 适 用 的 工 况 、 介 质 不 同 ， 超 临界 传 热 的 物理 
机 理 还 需要 深入 进行 实验 研究 。 
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